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We report an investigation of criticality in QCD at finite isovector chemical potential, fi^, and at 
zero temperature. At the critical point, « m^, we find that an uncharged scalar and pseudoscalar 
and a charged pseudoscalar meson become massless within the resolution of our measurement. The 
effective long distance theory therefore breaks 0(4) symmetry by charged pion condensation. This 
results in a rising quark number susceptibility. The baryon remains massive, as indicated by a 
vanishing baryon number susceptibility. 
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In the last few years there has been a resurgence of 
interest in the phase diagram of QCD 0. Perturbation 
theory and effective chiral models have been used to un- 
cover a rich phase structure at finite chemical potential 
and low temperatures [Q. This region of non- vanishing 
baryon chemical potential remains out of reach of present 
lattice computations, due to the Fermion sign problem 
H . However, the same techniques, namely effective chiral 
models and perturbation theory, have also been used to 
investigate the phase structure for non-vanishing isovec- 
tor chemical potential [0. This new direction in cou- 
pling space is amenable to lattice computations, which 
can therefore test the various methods applied to com- 
putations at finite baryon chemical potential. 

It is clear that the Dirac operator, M, is non-negative 
if MT = QMQ^ 1 for some operator Q. For two-flavour 
QCD with isovector chemical potential, 



M = 











and 



Q = 75 <H>t 2 , 



(1) 

the action is positive definite provided each flavour sat- 
isfies 75-hermiticity (D* = 75D75) and the quark masses 
are degenerate [^j . The phase diagram of this model was 
investigated by chiral perturbation theory and by per- 
turbative QCD and a second order phase transition in 
the 0(4) universality class was found at /13 = due to 
the condensation of charged pions pjj. A first quenched 
lattice investigation of this model has been reported re- 
cently Q and a phase transition due to pion condensation 
has indeed been found. When det M is real and positive, 
quenching does not produce the pathologies which can 
otherwise arise. 

We show here that it is possible to determine the sym- 
metry of the problem, and through universality then pre- 
dict the critical exponents, by investigating the spin- 
flavour content of the massless particles in the theory. 
Our quenched computations verify that a second order 
phase transition occurs at /J3 = m n . At, and above, the 
critical point the scalar tr, the neutral pion and one linear 
combination of charged pions becomes effectively mass- 



less. This lends support to the idea that the effective 
field theory should have an 0(4) symmetry broken down 
by giving mass to a charged pion, i.e., by pion conden- 
sation. If one is so inclined, this can be further tested 
by the usual, more expensive, method of finite-size scal- 
ing. However, as pointed out in this is non-trivial 
since higher dimensional operators in the lattice theory 
are important at the cutoffs which are presently accessi- 
ble. These could easily influence the effective exponents 
seen. Scaling studies become then many-fold more dif- 
ficult due to the necessity of taking finer lattices before 
investigating finite volume effects. 

The partition function of two-flavour QCD is 



VUe 



detM(m u , det M(m d ,fi d ), (2) 



where Sw is the Wilson gauge action and each determi- 
nant corresponds to one flavour of quarks. At the mo- 
ment we are interested in the case when m u = rrid = m. 
The chemical potentials for each flavour can be com- 
bined into the baryon and isovector chemical potentials 
fJ-Q = 3(fJ-u + fJ-d)/2 and fi 3 = (fi u - Hd)/2 respectively. It 
is instructive to invert these relations and write 



fJ-o 



fJ-3, 



Mo 



M3- 



(3) 



We shall work with /j,q = and /J3 > 0. Since the identi- 
fication of u and d quarks is arbitrary in the absence of 
electroweak interactions, the difference between positive 
and negative /U3 is purely conventional. A flip in the sign 
of /x 3 , as preferred in [Q, can be accomplished by just 
interchanging our definitions of u and d quark flavours. 
Quantum number densities are defined as 



T\ d\nZ 



(4) 



for i = and 3. These are linear combinations of the 
quark number densities. The quantity 713 is an order 
parameter for pion condensation. Other order parame- 
ters are the quark condensates (ipip) and (^75^2^). At 
the expected phase transition (tpi/i) would depart from 



2 



its vacuum value and (ipj5T2ip) would become non-zero. 
The added advantage in using is that it also flags the 
value of /is at which saturation sets in, with every site 
of the lattice being filled to capacity with quarks, and 
continuum physics can no longer be extracted. 

The quark number susceptibilities Q, Xij — drii/d/ij 
can also be used to probe the transition. The flavour 
off-diagonal susceptibility is given by 



T 



txM' u M-HvM' d M^ x ) 



- (trA^M- 1 ) {trM' d M- 1 ) 



(5) 



This differs from the usual fif = expression || through 
the subtraction of the last term, corresponding to non- 
vanishing quark densities. When a charged pion conden- 
sate forms, Xud must become non-zero and negative. Ex- 
pressions for the flavour diagonal susceptibilities can be 
read off from Q and modified by appropriate subtraction 
of the non-vanishing quark number densities. As long as 
the baryon remains massive, xo — dno/dfio should re- 
main zero. X3 is a combination of xo and Xud and hence 
gives no extra information. 

We explicitly use two staggered quark fields — one for 
each component of the Fermion matrix in eq. (Q). Wc 
adopt and extend the usual "cure" for Fermion doubling, 
and write the two-flavour (one up and one down) deter- 
minant as 



det M = (det M u M d ) 



1/4 



(6) 



Consequently, each trace in eqs. (fi]J|) comes with a factor 
of 1/4 multiplying it. There are two ways in which quark 
masses enter the problem — the sea quark mass appears 
in the determinants in eq. (|^) and influence the path 
integral measure, while the valence quark mass enters 
into the operators such as in eq. (||). In the quenched 
approximation the determinants are set to unity, so the 
sea quarks drop out of the problem but not the valence 
quarks. 

It might seem that the use of staggered quarks is a 
needless complication since the number of mesons of any 
given spin is much larger than that expected in two 
flavour QCD. However, only four of these are Goldstone 
pions at T = \i = 0. The others are split from them at 
any finite lattice spacing. This trick gives us the ability to 
bypass some of the problems seen at coarse lattice spac- 
ing. On finer lattices the number of Goldstone modes 
for staggered quarks is in any case a problem in defin- 
ing two staggered flavours and some other formulation of 
lattice quarks is probably better suited to the problem. 
One other advantage to using two independent staggered 
quark fields is that it allows us to give them unequal 
masses if we so wish. 

The partition function of eq. (Q) can be expressed 
through a transfer matrix — Z = TrT N (T, /zq, /K3). The 



transfer matrix T propagates field information from one 
slice of the lattice to another, and N is the number of 
slices. There may be several inequivalent ways to slice the 
lattice, giving rise to possibly inequivalent transfer matri- 
ces, but the partition function is unique. The eigenvalues 
of T, Xo > Ai > • • •, define masses (or screening masses) 
through the formula 1/m^ = log(Ao/Ai). A critical point 
implies massless modes for all the transfer matrices. We 
utilise the transfer matrix in the time direction and deter- 
mine the spin-flavour content of massless mesons which 
propagate in the time direction. 

Since the theory at finite chemical potential has a con- 
jugate quark ||, the symmetries of the transfer matrices 
must take this into account. In the problem at hand, 
there are two staggered quark fields which are conjugate 
to each other. The symmetry group at T = fi = there- 
fore contains an extra SU{2) factor corresponding to ar- 
bitrary rotations between these two fields. This symme- 
try is broken at finite [i. Recall that a chemical potential 
enters the lattice Dirac operator, M, through the time 
component of the covariant derivative because every tem- 
poral gauge field is multiplied by the fugacity, exp(a/i). 
As a result, for finite /X3, the Z{2) group generated by 
t — > — f does not commute with the new SU(2). The 
symmetry of the transfer matrix can be constructed by 
flipping the roles of the two quarks when flipping the 
sign of time. The irreps can be constructed by taking Z2 
invariant combinations such as 7r^. = uj^u ± djs d and 
7r~ = u"/sd or 7r + = dj^u. When electroweak corrections 
are neglected, as is the case here, one can also use the 
linear combinations n± — uj^d ± dj^u. The irrep labels 
of p^] can be retained with a reinterpretation of the old 
time-reversal quantum number as the product of time- 
reversal and the Zi flip of the up and down quarks above 
llTI. 




FIG. 1: The correlators for the dd part of 7r° (circles) and a 
(boxes) for a/X3 = 1 on a 24 x 8 3 lattice. The lines are fits to 
the spin projected part of the uu projection and rotated in 
flavour according to the symmetry of the transfer matrix. 

The importance of identifying the irreps of the trans- 



3 



fer matrix follows from an observation made in |l3|]: on 
any finite lattice an unphysical mode of propagation for 
a meson is that a quark and an antiquark travel in oppo- 
site directions around the lattice. This finite size effect 
is rendered innocuous at /i = by proper importance 
sampling. However, a chemical potential exacerbates 
this problem by encouraging forward propagation of the 
quark and backward propagation of the same flavour of 
the anti-quark. The correlation function of a pure flavour 
state is no longer symmetric in t, but has the form 



C(t) = Ae~ mt + A'e-™^"*-*), 



(7) 



with m 7^ m', which seems to indicate the absence of a 
transfer matrix. However, by recombining flavours into 
the irreps of the symmetries of the transfer matrix we 
recover the correct symmetric form. 

In Figure |l| we show an example of this. Note that 
masses fitted to such pure flavour correlators by the for- 
mula in eq. (0) are still useful, since the transfer matrix 
eigenstates (which are linear combinations of flavours) 
have the same masses. Tests of mass fits therefore include 
the check that the roles of the parameter pairs (A, m) and 
[A! ,ml) are flipped when fitting the uu and dd correla- 
tors, and that the eigenstates of the transfer matrix do, 
in fact, give the lower of these two masses. Figure |] does, 
in fact, show one such test, since the correlation functions 
plotted are the dd type, and the fits shown are plotted 
using parameters obtained from the uu type and flipped. 
We have further checked that local masses also show a 
plateau close to this mass. 

We have simulated quenched QCD on 8 4 and 12 4 lat- 
tices at j3 = 5.8941, which is the critical coupling for 
the finite temperature transition in quenched QCD with 
6 time slices. The lattice spacing is then fixed to be 
a = 1/6T C . In the quenched theory, since T c / Aj^ = 
1.15±0.05 @, the lattice cutoff is a" 1 = (6.9±0.3)A W 
The two lattice sizes differ by 50% in linear extent. 
This is necessary to see and control gross finite volume 
distortions, but not sufficient to perform a study with 
the finesse required to extract critical indices from fi- 
nite size effects. We have generated 40 stored config- 
urations for analysis on the smaller lattice and 70 on 
the larger. These are obtained with a Cabbibo-Marinari 
pseudo-heatbath algorithm where each SU(2) subgroup 
imbedded in SU (3) is touched thrice during every link 
update. An initial 5000 sweeps are discarded for ther- 
malisation and each subsequent stored configuration is 
separated from the previous one by 1000 sweeps. The 
successive stored configurations are therefore fully decor- 
related from each other, rendering the analysis of mea- 
surement errors fairly simple. 

Our methods for measuring the number densities and 
susceptibilities are the same as in ||. The chiral con- 
densate is measured from the same matrix inversions. 
Correlators are obtained with a point source and pro- 
jected on to zero 3-momentum by summing over the sink. 



The spectrum measurements use standard techniques for 
staggered quarks, including fitting to and extracting lo- 
cal masses from spin-projected parts of the correlators. 
Since we work on finite lattices, we have a finite reso- 
lution for the identification of massless modes. When a 
mass reaches the value l/aN t , the corresponding parti- 
cle becomes effectively massless and contributes to the 
long-distance effective theory. 

Since this is a first survey, we chose to work with a 
reasonably heavy valence quark mass — ma =1/6. At 
T = jti3 = this gave m^a — 1.065 ± 0.003 and m p a — 
1.323 ± 0.007. The baryon mass was m p a = 1.59 ± 0.05. 
These measurements are compatible with previous com- 
putations at nearby couplings |l5| . Since m^L « 12 on 
our larger lattice, we expect finite volume distortions to 
be small at this quark mass except in the neighbourhood 
of any critical point, where, in any case, we expect to see 
other physics as already explained. 
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FIG. 2: The order parameters. Boxes denote measurements 
on 8 4 lattices and circles those on 12 4 lattices. 



Measurements of the order parameters as functions of 
/13 are shown in Figure 0. We draw attention to two 
salient facts. First, (ifri/j) and both depart from their 
vacuum value at the same a\x\ w 1 k am^, and, second, 
the two lattice sizes yield identical results for the critical 
From the flattening of 713 at 0,^3 « 2 it is clear that 
saturation sets in and continuum physics can no longer 
be extracted. The critical point is reached far before 
this. The susceptibility \ud is zero at T = ^3 = [||. It 
begins to depart from this value, becoming negative, soon 
after fj,^. Its behaviour is closely related to the massless 
combination of charged pions. Throughout this range of 
/i3, no is consistent with zero, indicating that the baryon 
remains massive [ fl6|| . 

Figure ^ shows the /^-dependence of masses of the 7r° 
and 0+ obtained from zero-momentum correlators. The 
masses of the + and 7r° drop rapidly from their T = 
values and become nearly massless close to a/Zg. They 
remain massless till saturation sets in. The 7r+ mass is 
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FIG. 3: The spectrum of 7r+ (circles), cr+ (boxes) and 7r+ (di- 
amonds) obtained from temporal correlators on 12 4 lattices. 
The horizontal dotted line indicates the mass resolution avail- 
able on a lattice of this size. 



the most interesting. It is non-zero at [13 = due to the 
finite quark mass and degenerate with the ir+ . It does not 
change appreciably till /i§ , at which point it drops rapidly 
and becomes effectively massless soon after. The higher 
of the two masses in the flavour projected charged pion 
sector, corresponds to 7r_. This remains non-zero in the 
whole range of \i studied here. All the critical behaviour 
seen in the order parameters is driven by this. The rise 
in meson masses in the unphysical region a/is > 2 is 
easy to understand: propagation of mesons is hard on 
a saturated lattice. This is also observed in a free field 
theory computation at large /X3. 

A few caveats need to be stated. At the resolution 
available at these lattice sizes, a/i§ can be observed to 
lie between 1 and 1.25, as seen in Figure |[ Figure || 
shows that the mass resolution available on the 12 4 lat- 



pansions, and a lattice computation is the only possible 
tool to explore this region. The extraction of the phase 
diagram in these qualitatively different regions of ^3 is a 
major task, which we leave to the future. 

I would like to thank D. T. Son, M. Stephanov and R. 
Gavai for discussions. 
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of /z§ will be presented elsewhere. 
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